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1. Introduction 

The Bjorken-jc-dependence of parton distribution functions (PDFs) has been investigated by 
several lattice collaborations in recent years. In these studies, the lattice operators probing the 
nucleon always represent a local continuum operator. Spatial separations between the quark an- 
nihilation and creation operator on the lattice only appear in the context of discretized covariant 
derivatives. 

In this study, we explore the possibility to examine the dependence on the intrinsic parton 
momentum kj as well. For the purpose of calculating transverse momentum dependent parton 
distribution functions (TMDPDFs), we investigate non-local operators, constructed from quark 
and antiquark fields which are spatially separated. Gauge invariance is ensured by introduction of 
a straight gauge hnk (Wilson line) connecting the quark fields. 



2. TMDPDFs in SIDIS experiments 



In experimental processes like semi-inclusive deep inelastic scattering (SIDIS), one factorizes 
dominant diagrams into hard, perturbative processes and soft, non-perturbative parts [|l|, ^, rep- 
resented as shaded areas in fig. |l|a). The lower soft part can be parameterized in terms of nucleon 
TMDPDFs and corresponds to a correlator 

<I>[rl(fc,/',5) ^ l^^e'" m (f{i) \P,S) , (2.1) 
where = q{0)r^oj] (2-2) 

Here \P,S) represents a nucleon state of momentum P and spin 5, F is a Dirac matrix and k is the 
quark momentum. '^[o,£] is a gauge link connecting the quark operators. For SIDIS, the gauge link 
runs to infinity and back, as illustrated in fig. [I]b), see, e.g., ref. [Q]. 




Figure 1: a) Factorized tree level diagram of semi-inclusive deep inelastic scattering (SIDIS), b) Gauge 
link to infinity and back as in SIDIS, c) Straight gauge link 
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3. First approach to TMDPDFs from the lattice 

The idea of this study is to extract TMDPDFs from matrix ele- 
ments {P,S\ \P,S) calculated directly on the lattice. So far we 
restrict ourselves to straight gauge link operators: in the continuum 
limit, the gauge Unk follows a straight path connecting quark annihi- 
lation and creation operators, see fig. |l|c). The matrix elements thus 
calculated do not correspond exactly to those relevant for SIDIS, but 
the resulting TMDPDFs have an interesting probabilistic interpretation 
[Q]. On the lattice, the gauge link '^[qj] in the operator ^[^^{i) is a product of link variables. Since 
we work in Euclidean space-time, we can only evaluate spatial quark separations £ = {0,i). There- 
fore we restrict our lattice operators to a single slice T in Euclidean time. For quark separation / 
which do not lie on the x,y or z. axis, we can approximate a straight connection with a step-Uke link 
path, see fig. ^ (inset). 

To extract a signal for nucleon matrix elements of the operators in eq. (2^), we form ratios of 
nucleon three-point and two-point functions 



Figure 2: Step-like link 
path for r= (6,3,0). 
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(3.1) 



where F^p' and F^p' are suitable nucleon spin projection matrices^, and where B{t,P) h a nucleon 
interpolating operator, composed of u and d quark operators. Nucleon source and sink are placed at 
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Figure 3: a) Evaluation of the three-point function in the numerator of eq. (3.1 ) on the lattice (schematic), 
here for an operator with d-quarks. Only one of the two possible connected contractions of quark fields 
is shown. All-to-all propagators are avoided by combining three of the quark propagators into a sequential 
propagator (dark area), b) Overview of quark separations £ in the x-y-plane used in this investigation. We 
have calculated all quark separations which lie on the x- or y-axis, up to a length of 20 lattice units. In the 
first quadrant, we have included all quark separations up to a length of 8 lattice units (inner grey circle) and 
a selection of longer ones. 



'LHPC usesr^P' = 5(1 + 74) andr^P' = ^{1 + Y4){1 + ir^n) 
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times fsource = and t^i^k, respectively. The transfer matrix formalism reveals that at a sufficiently 
large distance from source and sink (0 ^ T ^ ?sink)> the ratio Rr{T;P,£) becomes T-independent. 
This plateau value is directly related to the value of the matrix element {P,S\ ^iat(^) l^)-^)- 

For the evaluation of eq. (3J), we apply the standard technique as described, e.g., in ref. |^], 
based on products of propagators and sequential propagators as illustrated in fig. ||a). In the fol- 
lowing, we present results for isovector operators (q = u — d), where disconnected contributions 
are absent. 



4. Test setup 



In our first numerical test calculations, we use 84 MILC gauge configurations from the NERSC 
archive [||,^. The configurations were produced with the AsqTad improved staggered quark action 
with 2+1 flavors. The lattice dimensions are x T = 20^ x 64, with a lattice spacing a ^ 0.124fm. 
The quark masses are aniuj = 0.030, and am^ = 0.050. The gauge configurations have been HYP 
smeared and bisected in the temporal direction, and we have selected only the time slices ... 31. 
We are using unsmeared propagators and sink-smeared sequential propagators previously calcu- 
lated by the LHPC group for these chopped gauge configurations (see, e.g. [^). The propagators 
have been calculated using domain wall fermions, with the quark mass tuned such that the pion 
mass « 596 MeV is approximately equal to the Goldstone pion mass for the staggered sea 
quark action. The source-sink separation is fsink — ^source = 10. LHPC has calculated sequential 
propagators for two nucleon momenta P = (0,0,0) and P = (—1,0,0). The latter corresponds to a 
momentum of 500 MeV in physical units. 

We have explored a number of link paths in all directions. Fig. ||b) illustrates our selection of 
quark separations £ in the x-y-p\ane. For our test runs, we have chosen two Dirac structures in the 
operator: the vector case F = 74 and the axial vector case F = 7375. 
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Figure 4: Sample plateau plots: Ry{t;P,£) is plotted versus t. The horizontal line and the error band 
indicate the plateau value and its error, extracted from the points (marked red), at T = 4,5 and 6. a) Real 
part of Rr{T;P, £) for F = 74, nucleon momentum P = (0, 0, 0), and a link path five units long in x direction, 
i.e. |£|=5. b) Imaginary part of 7?r(T,P,^) for F = 7375, nucleon momentum 7^= (0,0,0), andtheUnk 
path shown in fig. ^ i.e. \£\ = 6.7. 
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5. Preliminary Results 

In fig. ^ we show two sample plots of the ratio Ry{t;P,£) versus T. We obtain clean plateaus. 
In order to extract the plateau value Rr{P, tj, we average over time slices at T = 4, 5 and 6. 

Fig. ^a) shows all results Ry{P, tj for F = 74 and P = (0, 0, 0) for the 263 evaluated Unk paths. 
The signal is quite good, even for longer quark separations. We find that the correlator primarily 
depends on the separation \l\ between quark annihilation and creation operator. 

In figures § b) (vector case, unpolarized) and ^ a) (axial vector case, polarized) we select only 
link paths lying in the jc-y-plane. Furthermore, we identify groups of link paths which transform 
into one another under rotation or reflection in the x-y-^\me, and take the group average. In fig. 
|6|b) we show an example at non-zero nucleon momentum. For quark separations 1 on the positive 
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Figure 5: a) Results for F = 74, P = (0, 0, 0). We plot Re/?y^ (P, t) for all link paths versus the separation |£| 
of quark creation and annihilation operator, b) Re/?y^ {P ~ 0, t) versus \t\ for link paths in the x-y-plane. 
Results for link paths which transform into one another under rotation or reflection have been averaged. 



Dashed turquoise curve: fit to the data with a single Gaussian function //i (|£|), see eq. (5.1). Solid red curve: 
fit with the superposition of two Gaussian functions H2{\1\). The parameters determined from the fits are 
listed in tables [l| and |[ 
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Figure 6: a) Im/?r versus \l\ for F = 7375 (axial vector), momentum P — (0,0,0) and Unk paths in the 
x-y-plane. Results for link paths which transform into one another under rotation or reflection have been av- 
eraged, b) ^Re {^r(^) +^r(— ^)} for F = 74, £ on the positive .jc-axis, and non-zero nucleon momentum 




f = (—1,0,0). In both plots, the data have been fitted with a single Gaussian function //i (|^|), see eq. (5.1) 
Parameters determined from the fits are listed in table [l|. 
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jc-axis, nucleon momentum P = (—1,0,0) and F = 74, we plot ^Re |/?r(^) +^r(— ^)|- Here we 
refrain from averaging over the whole x-y-plane, because the parameterization indicates that the 
value is not invariant with respect to the Unk direction. 

We have tested the following fit functions to parameterize the \£\ dependence: 

HM :=Ciexp , //2(KI) :=Ciexp +C2exp . (5.1) 

The resulting fit parameters are listed in tables |I] and ^. We observe that the data can be well 
described by the sum of Gaussians in H2. 



6. A first glimpse of TMDPDFs from the lattice 

We can use the parameters of the fits in fig. |5]b) to calculate the first x-moment {n = 1) of the 
intrinsic transverse momentum dependence of the parton density f\ : 

fr'Ckr) ^ j\dx fjk- (I 2^ (^"^1 ^(0)7+%£]^(^) \P,S)^ ■ (6.1) 
It turns out that 

f&kr) = J ^e-'^'r-^-Ry,{P = 0, \It\) , (6.2) 
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Table 1: Fit parameters determined from the single Gaussian fits with Hi shown in fig. |^b) and || 
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Table 2: Fit parameters determined from the double Gaussian fit with H2 shown in fig. ^b) 
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Figure 7: Upper curve: /"j^J (^r), calculated from a Fourier transform of Ry^{P = 0, |^|). The upper curve 
is the sum of the two dashed curves, which show the Fourier transforms of the two Gaussian components in 
fit function H2{\t\)- The results are not renormaUzed. 
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where "lat" indicates that the operator has not been renormaUzed. For the single Gaussian fit (func- 
tion H\), we find a root mean square transverse momentum of \J (kj) = 2/o = (563 ± 12)MeV, 
see table |. This is very well compatible with a value of k, 500 MeV used in recent phe- 
nomenological investigations of HERMES data on SIDIS [|^], based on the factorized Ansatz 
fi{x,kT) = fi (■^)exp[— ^^/(^|.)]/[7r(^j)]. Note however, that such a comparison with phenomeno- 
logical results has to be taken with due caution, since the effect of renormalization of the non-local 
operators ^j^j could, in principle, affect the ^dependence of Ry- The result for the double Gaussian 
fit (function H2) is shown in fig. ^ and table ^ We obtain a root mean square of the transverse 

momentum of ^J^) = (702 ± 12) MeV. 

7. Conclusions and outlook 

We have calculated nucleon matrix elements {P,S\ g(0)r^[o.f] ^(^) 1^)'^) with a finite sepai^a- 
tion I of the quark operators. It turns out that the dependence on i is approximately Gaussian in 
the channels we explored. We have used our data to obtain a first, preliminary result on transverse 
momentum dependent parton distribution functions (TMDPDFs). The root mean square transverse 

momentum \J (kj) of our unrenormalized result for /"^^(kr) is (563 ± 12) MeV for a single Gaus- 
sian fit, a value which is compatible with phenomenological results. It will be interesting to study 
the correlations with respect to the quark separation £ and the nucleon momentum P. We also plan 
to investigate whether there is a lattice analogy to link paths extending to infinity and back. This 
would enable us to calculate the TMDPDFs directly relevant to phenomenology. 
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